Response to sub-threshold stimulus is enhanced by spatially heterogeneous activity

Sub-threshold stimuli cannot initiate excitations in active media, but surprisingly as we show in this paper, they can alter the time-evolution of spatially heterogeneous activity by modifying the recovery dynamics. This results in significant reduction of waveback velocity which may lead to spatial coherence, terminating all activity in the medium including spatiotemporal chaos. We analytically derive model-independent conditions for which such behavior can be observed.Comment: 5 pages, 5 figure

Selfsimilar solutions in a sector for a quasilinear parabolic equation

We study a two-point free boundary problem in a sector for a quasilinear parabolic equation. The boundary conditions are assumed to be spatially and temporally "self-similar" in a special way. We prove the existence, uniqueness and asymptotic stability of an expanding solution which is self-similar at discrete times. We also study the existence and uniqueness of a shrinking solution which is self-similar at discrete times.Comment: 23 page

Incorporating Inductances in Tissue-Scale Models of Cardiac Electrophysiology

In standard models of cardiac electrophysiology, including the bidomain and monodomain models, local perturbations can propagate at infinite speed. We address this unrealistic property by developing a hyperbolic bidomain model that is based on a generalization of Ohm's law with a Cattaneo-type model for the fluxes. Further, we obtain a hyperbolic monodomain model in the case that the intracellular and extracellular conductivity tensors have the same anisotropy ratio. In one spatial dimension, the hyperbolic monodomain model is equivalent to a cable model that includes axial inductances, and the relaxation times of the Cattaneo fluxes are strictly related to these inductances. A purely linear analysis shows that the inductances are negligible, but models of cardiac electrophysiology are highly nonlinear, and linear predictions may not capture the fully nonlinear dynamics. In fact, contrary to the linear analysis, we show that for simple nonlinear ionic models, an increase in conduction velocity is obtained for small and moderate values of the relaxation time. A similar behavior is also demonstrated with biophysically detailed ionic models. Using the Fenton-Karma model along with a low-order finite element spatial discretization, we numerically analyze differences between the standard monodomain model and the hyperbolic monodomain model. In a simple benchmark test, we show that the propagation of the action potential is strongly influenced by the alignment of the fibers with respect to the mesh in both the parabolic and hyperbolic models when using relatively coarse spatial discretizations. Accurate predictions of the conduction velocity require computational mesh spacings on the order of a single cardiac cell. We also compare the two formulations in the case of spiral break up and atrial fibrillation in an anatomically detailed model of the left atrium, and [...].Comment: 20 pages, 12 figure

Identification of IKr Kinetics and Drug Binding in Native Myocytes

Determining the effect of a compound on IKr is a standard screen for drug safety. Often the effect is described using a single IC50 value, which is unable to capture complex effects of a drug. Using verapamil as an example, we present a method for using recordings from native myocytes at several drug doses along with qualitative features of IKr from published studies of HERG current to estimate parameters in a mathematical model of the drug effect on IKr. IKr was recorded from canine left ventricular myocytes using ruptured patch techniques. A voltage command protocol was used to record tail currents at voltages from −70 to −20 mV, following activating pulses over a wide range of voltages and pulse durations. Model equations were taken from a published IKr Markov model and the drug was modeled as binding to the open state. Parameters were estimated using a combined global and local optimization algorithm based on collected data with two additional constraints on IKrI–V relation and IKr inactivation. The method produced models that quantitatively reproduce both the control IKr kinetics and dose dependent changes in the current. In addition, the model exhibited use and rate dependence. The results suggest that: (1) the technique proposed here has the practical potential to develop data-driven models that quantitatively reproduce channel behavior in native myocytes; (2) the method can capture important drug effects that cannot be reproduced by the IC50 method. Although the method was developed for IKr, the same strategy can be applied to other ion channels, once appropriate channel-specific voltage protocols and qualitative features are identified

Simulation and Mechanistic Investigation of the Arrhythmogenic Role of the Late Sodium Current in Human Heart Failure

Heart failure constitutes a major public health problem worldwide. The electrophysiological remodeling of failing hearts sets the stage for malignant arrhythmias, in which the role of the late Na+ current (INaL) is relevant and is currently under investigation. In this study we examined the role of INaL in the electrophysiological phenotype of ventricular myocytes, and its proarrhythmic effects in the failing heart. A model for cellular heart failure was proposed using a modified version of Grandi et al. model for human ventricular action potential that incorporates the formulation of INaL. A sensitivity analysis of the model was performed and simulations of the pathological electrical activity of the cell were conducted. The proposed model for the human INaL and the electrophysiological remodeling of myocytes from failing hearts accurately reproduce experimental observations. The sensitivity analysis of the modulation of electrophysiological parameters of myocytes from failing hearts due to ion channels remodeling, revealed a role for INaL in the prolongation of action potential duration (APD), triangulation of the shape of the AP, and changes in Ca2+ transient. A mechanistic investigation of intracellular Na+ accumulation and APD shortening with increasing frequency of stimulation of failing myocytes revealed a role for the Na+/K+ pump, the Na+/Ca2+ exchanger and INaL. The results of the simulations also showed that in failing myocytes, the enhancement of INaL increased the reverse rate-dependent APD prolongation and the probability of initiating early afterdepolarizations. The electrophysiological remodeling of failing hearts and especially the enhancement of the INaL prolong APD and alter Ca2+ transient facilitating the development of early afterdepolarizations. An enhanced INaL appears to be an important contributor to the electrophysiological phenotype and to the dysregulation of [Ca2+]i homeostasis of failing myocytes

Immune cell counts and risks of respiratory infections among infants exposed pre- and postnatally to organochlorine compounds: a prospective study

<p>Abstract</p> <p>Background</p> <p>Early-life chemical exposure may influence immune system development, subsequently affecting child health. We investigated immunomodulatory potentials of polychlorinated biphenyls (PCBs) and <it>p,p'</it>-DDE in infants.</p> <p>Methods</p> <p>Prenatal exposure to PCBs and <it>p,p'</it>-DDE was estimated from maternal serum concentrations during pregnancy. Postnatal exposure was calculated from concentrations of the compounds in mother's milk, total number of nursing days, and percentage of full nursing each week during the 3 month nursing period. Number and types of infections among infants were registered by the mothers (N = 190). White blood cell counts (N = 86) and lymphocyte subsets (N = 52) were analyzed in a subgroup of infants at 3 months of age.</p> <p>Results</p> <p>Infants with the highest prenatal exposure to PCB congeners CB-28, CB-52 and CB-101 had an increased risk of respiratory infection during the study period. In contrast, the infection odds ratios (ORs) were highest among infants with the lowest prenatal mono-<it>ortho </it>PCB (CB-105, CB-118, CB-156, CB-167) and di-<it>ortho </it>PCB (CB-138, CB-153, CB-180) exposure, and postnatal mono- and di-<it>ortho </it>PCB, and <it>p,p'</it>-DDE exposure. Similar results were found for pre- and postnatal CB-153 exposure, a good marker for total PCB exposure. Altogether, a negative relationship was indicated between infections and total organochlorine compound exposure during the whole pre- and postnatal period. Prenatal exposure to CB-28, CB-52 and CB-101 was positively associated with numbers of lymphocytes and monocytes in infants 3 months after delivery. Prenatal exposure to <it>p,p'</it>-DDE was negatively associated with the percentage of eosinophils. No significant associations were found between PCB and <it>p,p'</it>-DDE exposure and numbers/percentages of lymphocyte subsets, after adjustment for potential confounders.</p> <p>Conclusion</p> <p>This hypothesis generating study suggests that background exposure to PCBs and <it>p,p'</it>-DDE early in life modulate immune system development. Strong correlations between mono- and di-<it>ortho </it>PCBs, and <it>p,p'</it>-DDE exposures make it difficult to identify the most important contributor to the suggested immunomodulation, and to separate effects due to pre- and postnatal exposure. The suggested PCB and <it>p,p'</it>-DDE modulation of infection risks may have consequences for the health development during childhood, since respiratory infections early in life may be risk factors for asthma and middle ear infections.</p

Effects of Electrical and Structural Remodeling on Atrial Fibrillation Maintenance: A Simulation Study

Atrial fibrillation, a common cardiac arrhythmia, often progresses unfavourably: in patients with long-term atrial fibrillation, fibrillatory episodes are typically of increased duration and frequency of occurrence relative to healthy controls. This is due to electrical, structural, and contractile remodeling processes. We investigated mechanisms of how electrical and structural remodeling contribute to perpetuation of simulated atrial fibrillation, using a mathematical model of the human atrial action potential incorporated into an anatomically realistic three-dimensional structural model of the human atria. Electrical and structural remodeling both shortened the atrial wavelength - electrical remodeling primarily through a decrease in action potential duration, while structural remodeling primarily slowed conduction. The decrease in wavelength correlates with an increase in the average duration of atrial fibrillation/flutter episodes. The dependence of reentry duration on wavelength was the same for electrical vs. structural remodeling. However, the dynamics during atrial reentry varied between electrical, structural, and combined electrical and structural remodeling in several ways, including: (i) with structural remodeling there were more occurrences of fragmented wavefronts and hence more filaments than during electrical remodeling; (ii) dominant waves anchored around different anatomical obstacles in electrical vs. structural remodeling; (iii) dominant waves were often not anchored in combined electrical and structural remodeling. We conclude that, in simulated atrial fibrillation, the wavelength dependence of reentry duration is similar for electrical and structural remodeling, despite major differences in overall dynamics, including maximal number of filaments, wave fragmentation, restitution properties, and whether dominant waves are anchored to anatomical obstacles or spiralling freely